Emergence of weight-topology correlations in complex scale-free networks 
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Different weighted scale-free networks show weights-topology correlations indicated by the non 
linear scaling of the node strength with its connectivity. In this paper we show that networks with 
and without weight-topology correlations can emerge from the same simple growth dynamics of 
the node connectivities and of the link weights. A weighted fitness network is introduced in which 
both nodes and links are assigned intrinsic fitness. This model can show a local dependence of the 
weight-topology correlations and can undergo a phase transition to a state in which the network is 
dominated by few links which acquire a finite fraction of the total weight of the network. 
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Recently many complex systems have been described 
in terms of the underline network structure responsible 
for making one global entity performing specific functions 
from many otherwise independent elements 00>l5]- The 
scale-free degree distribution is a relevant characteristic 
of a number of complex networks which originates, for 
a large number of real networks, from their growing na- 
ture and in the "preferential attachment" of new links 
to well connected nodes. The minimal model which ac- 
counts for these properties is the Barabasi-Albert (BA) 
model Nevertheless since the early days of research 
in complex scale-free networks it has been recognized 
that the links are not equivalent like in the BA model 
but they have a different weight Wij for every link 
For example in the Internet 0, @ every connection has 
a given bandwidth, in the coauthorship network ev- 
ery couple of coauthors has a different number of pa- 
pers together, in the airport network every connec- 
tion between two airports has a different traffic and in 
the shareholder network the investments involve dif- 
ferent volumes of shares. An important quantity to char- 
acterize the relevance of a node in a weighted network is 
its strength defined as the sum of all the weights of 
incoming and outgoing links Sj = Ylj w i,j- Two classes 
of weighted networks can be identified looking at the be- 
havior of the strength of the node with the connectivity 
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with 8 > 1. In real systems an example of class I net- 
works is the coauthorship network0 while example of 
class II networks are the airport networks the share- 
holder networks and the Internet 0. The proposed 
classification discriminate between networks in which the 
connectivity of the node doesn't affect the weights of 
its links (class I networks) and networks in which the 
connectivity strongly influence them (class II networks). 
While studying specific systems one could find charac- 
teristic feature of them, a common behavior of different 
networks indicate an universality which ask for an expla- 
nation. The present models of weighted networks in the 



literature generate class I networks 0, El E3 and class II 
networks 12, 

m by different mecha^s ln particular 
the non linear scaling of strength and connectivity is as- 
sumed as an hypothesis in Ref. ^3] and it is the signature 
of some non linear coupling between weight and strength 
in Ref. 01 . In this paper we want to show that it is 
possible to explain the emergence of class I and class II 
networks under a common framework. We first propose 
a very schematic model, the preferential strengthening 
growing network model. In this model both the weights 
of the links and the connectivities of the nodes grow in 
time: networks of type II emerge only when the rate m' at 
which weights are strengthened is higher than the rate m 
at which new links are established. Then we propose the 
weighted fitness model where a fitness variable is assigned 
to each node and to each link of the network accounting 
for the intrinsic capacity of nodes to acquire new links 
and of links to acquire more weight. This model gener- 
ates networks of class I or class II depending on the value 
of the ratio m'/m. Moreover in the fitness weighted net- 
work we can distinguish two phases as a function of the 
ratio m'/m. In the first phase all links have an infinites- 
imal amount of the total weight of the network while in 
the second phase few connections grab a finite fraction 
of all the weight in the network. 

The basic idea of both models is that in many sys- 
tems the weight of the links is evolving in time while 
the network grows. For example in the coauthorship net- 
work two scientists in time increase the number of papers 
written together, in the airport network, while more con- 
nections are established between new airports and exist- 
ing ones, the traffic increases for all the existing connec- 
tions, in the shareholding network the volume of invest- 
ment change with time. Consequently in our model as 
the network grows with new nodes joining the network, 
also the total weight of the network increases and the 
weight of the existing connections is strengthened. For 
simplicity we assume that the two processes appear at 
the same time scale but with different rate. We have as- 
sumed a "preferential attachment" [4j mechanism for the 
strengthening of the connections of the nodes, i.e. links 
will increase their weight in proportion of their actual 
weight. 
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The preferential strengthening growing network - Given 
a starting network we assume that at each time a new 
node arrives in the network and it attaches to links to 
the existing nodes of the network establishing connec- 
tions with a minimal weight wq ■ At the same time other 
ml links increase their weight of an equal amount Wq. 
These links are chosen following a preferential attach- 
ment driven by the weight of the links. In particular at 
each time step i) we add a node and to links following 
the BA model 4]; it) we choose to' links and increase 
their weight by wq. We perform the choice of the links 
first by choosing a node with probability proportional to 
its strength IT = Si/J2j s j> an d then by choosing one 
of its links with probability proportional to its weight 
. The mean field p3| dynamic equa- 



lly = wij/E^ 

tion for the weight of a link is: 
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The solution of this equation is Wi.j = wq {t/tij) a where 
fj.j = max(ij,tj) and ti is the time at which 



node i was added to the network. The dynamics of the 
connectivities is the simple BA |4j model with ki(t) fol- 
lowing, in mean field, dki(t)/dt = rnki(t)/Y2j kj(t) hav- 
ing as solution ki(t) = m(t/ti) 1//2 . If we define ptj 
the probability that a node i is connected to node j 
(pij — m/(2y/titj), the strength of a node Si is given by 



Si(t) 



Asymptotically in time we found 



ki if a < 1/2 ml < to 

.s; v < fe,Tog(fe,) if a = 1/2 m' = to 
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if a > 1/2 ml > to 



(2) 



with the the weights satisfying Wij(t) oc (k m i n ) 2a where 
k m in = min(fc,', kj) making of this model when to' < to 
a class I network and when to' > to a class II network. 
Consequently just with this simple model we found that 
class II networks appear spontaneously when the rate of 
the strengthening of the weights is higher than the rate 
of addition of new links. Moreover it is easy to show 
P| that both the degree and the weight distributions are 

scale-free with P(k) = 2tm 2 k~ 3 and P{w) — tw^** j^Tt? 
where 6 = 1/a + 1. An important quantity to character- 
ize the inhomogeneities of the weights of the links end- 
ing on the same node is Y 2 = J2j(wij /si) 2 ^5[. When 
all the weights give a similar contribution Y 2 ~ 1/fej 
when there is one link with a weight Wi.j ~ s, one finds 
Y 2 ~ 0(\). In our preferential strengthening model we 
found relevant inhomogeneities also in the class I net- 
work for to' > to/3. In fact we have Y 2 (k) ~ \/k for 
to' < to/3 while Y 2 {k) ~ k 4a ~ 2 for to > to' > to/3. For 
class II networks ml > to we have Y 2 (k) ~ 0(1) . Wc 
have performed numerical simulations of the model that 
are well captured by the above mean-field treatment. In 
Fig^ we show the evolution of the strength of a node 
as a function of time for networks with a different ratio 
ml I'm. Moreover we show the non trivial dependence of 
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FIG. 1: Characteristic behavior of the basic model with 
a temporal power-law evolution of the strength Sj(t) (left 
graph) and power-law dependence of the strength of a node 
with connectivity (right graph). The data shown are for the 
preferential streghtening model with in ascending ordering 
(m,m') = (2,1), (2, 3), (1,3), N = 10 5 averaged over 100 re- 
alizations. 



s(k) = (si\ki — k) that deviates from linearity if to' > to 
as predicted by (J2J- 

The fitness weighted network- To justify how correla- 
tions between connectivities can emerge we consider the 
fitness weighted network. In this model a random fitness 
rji from a p(rj) distribution is assigned to each node of the 
network as in and the network grows with new links 
attached preferentially to high fitness and high connec- 
tivity nodes. Moreover a fitness £ij is assigned to each 
link (i, i) o f the network and as the network grows fol- 
lowing 16] links with higher fitness and higher weights 
are preferentially strengthened. The distribution from 
which is extracted the link fitness £ij will in general de- 
pend on the value of the fitness of the two linked nodes 
p(^i,j\T]i,rjj). For this fitness weighted model the weight 
of single links will increase following, in mean field, 
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Following the same technique used for the fitness model 
we assume that J2e v €t,t' w t,t' — 2m'woC't. With this 
ansatz the solution of this equation is straightforward 
and is given by 



w, 



Wq 



(4) 



with a(£ij) = ^jr, ti j = max(tj,tj*) where fcj indicates 
the time at with node i was added to the network. Ob- 
serve that C has the lower bond C < £m = max^j- 
since the weight of a link for construction cannot grow 
faster than linearly. For the dynamics of new connec- 
tions we use the fitness modelll6|. The evolution of the 
connectivity of a node i is then given by dki(t)/dt = 

mrjiki(t)/ Vjkj{t)j with solution k{(t) = m 

and C given by the self-consistent equation for cflrjf 
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1 = J dr]p(rj)rj/(C—T]). Substituting © in J2e,e> ii,i' w W 
and assuming that this sum is self-averaging, we get the 
equation for C" 



Tft 

— = I drjip(rn) I d€i tj p(£ itj \r}i) 
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Note that the right hand side of this expression has 
clearly its maximum for C = £m that is the minimum 
possible value of C. On the other side the strength of 
a node i, with fitness m can be well approximated by the 
value of its average on the fitness of the other nodes and 
the average on the fitness of its links. In this approxima- 
tion we have, asymptotically in time 



for ?7 . < ^ max . C .. 



(6) 



with 9{ri) > 1 depending on the distributions 
p( r ))jP{ii,j\ r liT r lj)- For every specification of the model 
one find the weighted scale-free network in class I or in 
class II depending on the value of m'/m. 

Let us for a moment consider the case of uncorrelated 
{f?},{£} distributions with p(rf) and p(£) uniform with 
i] G (0,t]m) and £ 6 (0,£m)- In this case we have that 
© holds with e(rf) = t M C/{C'ri). The average s(k) is 
then given by s(k) — J drjp(r))s v (k)P(k\rj) which gives 



(7) 



with = -C'/£m + C/r]M + 1 > 0. Since C is a mono- 
tonic function of the ratio r ~ m'/m this means that 
there will be a certain value tq such that for r < tq the 
network is in class I and for r > ro the network is in class 
II. As in the preferential strengthening network model 
here we found that Y 2 (k) indicates relevant heterogene- 
ity also in class I networks. In fact we have for class I net- 
works Y 2 (k) ~ 1/k for C > 2x M C/rj M and Y 2 (k) ~ fc~ A 
with A € (0,1) and A = min(C"/6/ + C/r] M ,-l - 
C"/(2£m) + C/m) for & E {ZmC/^mMm/vm)- For 
class II networks instead we found Y 2 (k) ~ O(l). The 
above treatment is valid as long as the solutions of © and 
the self-consistent equation for C exist, i.e. p(t]m) > 
and p(£m) > 0. When on the contrary p(rj) ~ (t]m — v) k 
with k > we have a phase transition (Bose-Einstein 
condensation on complex networks -BECCN)[ljjin which 
one node grabs a finite fraction of all the links. Similarly 
if p(£m) ~ (Cm — with k' > we have a "Bose- 
Einstein condensation" of the weights of the network, i.e. 
one link grab a finite fraction of all the weight of the net- 
work. The study of a simple case of unlimited support of 
the probability distributions p (v) = e~ v (or p(£) = e~^) 
can be easily be carried out [l(j-Here we observe that 
if both fitness of the nodes and fitness of the links are 
exponentially distributed, the scaling © remains valid. 

Even more interesting is the case in which the fitness of 
the links is functionally related to the fitness of the nodes 




FIG. 2: Power-law distribution for the strength in the 
weighted fitness model as a function of k for = r\i + r\j 
and uniform distribution of the r\ £ [0, 1] for the parameters 
(m,m') = (2, 1), (1, 10). The network have 10 nodes and the 
results are averaged over 100 realizations. Inset: Solution of 
© as a function of the ratio m'/m for uniform distribution of 
the r\ and fjj = r)i + r]j (solid line),£ij = r\ir\j (dashed line). 
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FIG. 3: Fraction of the total weight covered by the top 9 
links of the network in a weighted fitness model of m — 1 as 
a function of m! . The data are shown for a graph of N = 
10 3 nodes averaged over 500 realizations. The dashed line 
indicates the value m'/m = 3.28 over which IpjJ doesn't have 
a solution and we predict that few nodes grab a finite fraction 
of the total weight. 



£ij = f(r)i,r]j). This is a natural assumption in many 
real complex systems. For example, more productive 
scientists will produce more papers both with existing 
coauthors and with new collaborations, airports which 
attract more connections will reasonably also be the air- 
ports which attract more traffic, richer agents will invest 
more in new investments and in their present portfolio. 
Two relevant examples of shape of f(r]i, rjj) would be for 
example f(rji,r)j) = n + m and f(r]i,r}j) = rum. We will 
first consider the case of additive f(rji, m) = m + m with 
uniform distribution of the rj and r\ € (0,?7m)- In this 
case C/tim = 1.255 ... HU and C < 2rj M = xm ■ br Fig. 
Epnset) we plot the value of C as a function of the ratio 
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w! I'm as predicted by JSJ) in this case. In C = 2t]m the 
integral in Q5J converge allowing for the maximal value 
of the ratio ^ = r c = 3.28 ... for © to be valid. This 
scenario is indeed reminiscent of BECCN. In fact the 
convergence of © for C — 1t\m at m'/m = r c is the 
signature that for higher values of the ratio m'/m there 
will be some links that have a finite fraction of the overall 
weight of the network. We performed numerical simula- 
tions of the fitness weighted network model as a function 
of the ratio m'/m. In Fig. [3]we plot the fraction over the 
total weight of the weight of the top links of networks fo 
size N = 10 3 . We observe that for m'/m > r c — 3.28 . . . 
the fraction of weight acquired by few top nodes is a fi- 
nite fraction of the total weight. The difference between 
this transition and BECCN is that in the BECCN a sin- 
gle node grabs a finite fraction of the links while here 
we have few links that grab a finite fraction of the to- 
tal weight of the network. In fact we have a finite den- 
sity of links with maximal fitness p(£zj = 2t)m) > 0. 
Below this phase transition we found s(k) ~ k e with 

6 = 1 + (2C- C')/rjM for C > 2C,i.e. the network is a 
class II network for m'/m > ro = 1.18 . . . and s(k) ~ k 
for C < 2C, i.e. the network is a class I network for 
m'/m < ro. In Fig. [21 we report the average strength 
of node as a function of the connectivity k, s(k) for net- 
work in class I (m'/m — 0.5) and in class II (m'/m = 10). 
The degree and the weight distribution below the phase 
transition follows a power-law this time with a logarith- 
mic correction P(k) oc , k 7,. and P(w) oc -n^-. — ryr with 

V I log(fc) V I (log(™)) 2 

7 = 1 + C = 2.255 ... and 8 = 1 + C/2 while above the 
phase transition in the degree and weight distribution is 



show the presence of few condensed links. 

This scenario we are presenting is not restricted to 
the specific case £jj = r/i + r)j. For example the phase 
transition indicating the "condensation" of the weights is 
also present in the case of the multiplicative dependence 
of the fitness of the links with the fitness of the nodes 
£i j = ViVj arL d uniform p(rj) distribution with the criti- 
cal ratio r c ~ 1.334. In the Inset of Fig. |2Jwe report the 
dependence of C on the ratio m'/m for this last case. In 
this case the network is a class I network for C > Ct)m 
and a class II network for C < Cr\u with = Ctjm/C 
and the critical value of m'/m = tq = 0.59. 

In conclusion we presented a general framework for 
weighted networks which describe the emergence of 
weight-topology correlations. In this framework both 
node connections and link weights increase following a 
preferential attachment rule. The networks of class I or 
class II (with significant weight-topology correlations) are 
obtained as a function ratio m'/m between the rate m' at 
which links are strengthened and the rate m at which new 
connections are established. Moreover also in class I net- 
works there can be inhomogeneities in the weights of node 
edges reflected in a non trivial Y 2 (k). For the weighted 
fitness model which accounts for networks with degree- 
degree correlations the hetereogenieties in the weights of 
the nodes can be so strong that for m'/m > r c there is a 
structural phase transition to a phase in which few links 
grab a finite fraction of the total weight of the network. 
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